Microscopic Calculation of Pre-Compound Excitation Energies 
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We introduce a microscopic approach for calculating the excitation energies of systems formed 
during heavy-ion collisions. The method is based on time-dependent Hartree-Fock (TDHF) theory 
and allows the study of the excitation energy as a function of time or ion-ion separation distance. 
We discuss how this excitation energy is related to the estimate of the excitation energy using the 
reaction Q-value, as well as its implications for dinuclear pre-compound systems formed during 
heavy-ion collisions. 
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During a heavy-ion collision the initial kinetic energy 
of the two nuclei is gradually converted into internal ex- 
citation. This excitation may result in exciting various 
types of collective modes or intrinsic excitations that lead 
to the heating of the combined nuclear system. The 
excitation of projectile-like and target-like fragments in 
deep-inelastic heavy-ion collisions is a well known man- 
ifestation of this phenomenon, which has been studied 
theoretically and experimentally On the other hand, 
for a dinuclear pre-compound system formed during a 
heavy-ion collision, the mode of decay may critically de- 
pend on the excitation energy of the system. Examples 
include collisions that may be candidates for the forma- 
tion of superheavy elements in hot or cold fusion reac- 
tions 0, 01 . Excitation energy is also an important in- 
gredient for the stochastic mean-field approach to nu- 
clear dynamics, which deals with fluctuations of collec- 
tive motion in addition to the average evolution [1, Q. 
While exclusive measurements of excitation energy may 
be possible for equilibrated systems (e.g. compound nu- 
clei formed in complete fusion or fragments produced in 
deep-inelastic collisions), the intermediate states formed 
during a collision have a short life time and are not ex- 
pected to be fully equilibrated thus making the measure- 
ment as well as the interpretation very difficult Q. For 
all these reasons it is desirable to develop a dynamical 
approach for calculating the excitation energy of the sys- 
tem as it evolves in time. 

It is generally acknowledged that the TDHF theory 
provides a useful foundation for a fully microscopic many- 
body theory of low-energy heavy-ion reactions 0, Q- 
While the long-time evolution in TDHF theory may not 
be very reliable, recent three-dimensional TDHF calcula- 
tions with no symmetry assumptions and using modern 
Skyrme forces have shown to accurately reproduce phe- 
nomena determined by the initial stages of the heavy- 
ion dynamics [§, [l^, Recently we have developed 
the density-constrained TDHF (DC-TDHF) method [ll, 
which is based on the generalization of the density con- 
straint method developed earlier [l^. We have shown 



that using the DC-TDHF method ion-ion potential bar- 
riers can be accurately produced [13, [l^, as these 
calculations also depend on early stages of the ion-ion 
dynamics. In addition, one-body energy dissipation ex- 
tracted from TDHF for low-energy fusion reactions was 
found to be in agreement with the friction coefficients 
based on the linear response theory as well as those in 
models where the dissipation was specifically adjusted to 
describe experiments AH of these new results sug- 

gest that TDHF dynamics provide a good description of 
the early stages of heavy-ion collisions. 

In this manuscript we outline a microscopic approach 
for calculating excitation energies of composite or din- 
uclear systems formed during heavy-ion collisions. The 
goal of the approach is to divide the TDHF motion into 
a collective and intrinsic part. The major assumption in 
achieving this goal is to assume that the collective part is 
primarily determined by the density pijc, t) and the cur- 
rent j(r,t). Consequently, the excitation energy can be 
formally written as 



E*{t) = Etdhf ~ E,oU ip{t),i{t)) , 



(1) 



where Etdhf is the total energy of the dynamical sys- 
tem, which is a conserved quantity, and Ecou represents 
the collective energy of the system. In the next step we 
break up the collective energy into two parts 

Ecoii (t) = Eku, {p{t),m) + Edc {p{t)) , (2) 
where Ekin represents the kinetic part and is given by 
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(3) 



which is asymptotically equivalent to the kinetic energy 
of the relative motion, ^pR^, where p is the reduced 
mass and R{t) is the ion-ion separation distance. The 
energy Edc is the lowest-energy state of all possible 
TDHF states with the same density and is required to 
have zero excitation energy. This state is found by using 
the density-constraint method [l^ [l^ , which minimizes 
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the energy while holding the instantaneous TDHF den- 
sity constant. We have previously shown JJj that Ejjc 
equals the ion-ion potential, V{R), shifted by the bind- 
ing energies of the participating nuclei, which allows us 
to write 

Ecoiiit) = EkMt),m) + ViRit)) + Ea, + Ea, , (4) 

where Ea^ and Ea2 denote the Hartree-Fock energies 
calculated for the two nuclei using the same effective in- 
teraction. The dynamics of the ion-ion separation R{t) 
can be extracted from an unrestricted TDHF run. Using 
E*(t) and R{t), we can deduce the excitation energy as 
a function of the distance parameter, E*{R). 

The computation of the excitation energy as outlined 
above is numerically very intensive, primarily due to the 
density-constraint calculation. Calculations were done in 
3-D geometry and using the full Skyrmc force (SLy4) [l^ 
without the center-of-mass correction as described in 
Ref. The numerical accuracy of the static bind- 

ing energies and the deviation from the point Coulomb 
energy in the initial state of the collision dynamics is 
on the order of 50 — 200 keV. We have performed den- 
sity constraint calculations at every 20 fm/c. For the 
calculation of the ion-ion separation distance R we use 
the hybrid method, which relates the coordinate to the 
quadrupole moment for small R values, as described in 
Ref. jl6|]. The accuracy of the density constraint calcu- 
lations is commensurate with the accuracy of the static 
calculations. 
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FIG. 1: Internuclear potential V{R) for the head-on collision 
of the ^^O+^^O system for various -Ec.m. values. The relative 
ground-state binding energy of the "^^S nucleus is represented 
by the —Qgg value. 

In order to develop a better insight into the excitation 
energy given by Eq. ([1]) we have first studied two spheri- 
cal systems, ^^0-1-^^0 and '"'Ca-l-^^Ca. In literature one 
commonly defines the excitation energy for a particular 
reaction as 

E* — Ec.m. + Qgg , (5) 



where Qgg is the mass difference between the two ini- 
tial nuclei and the combined system in its ground state. 
While this expression is correct relative to the ground 
state of the composite system, it does not accurately de- 
scribe the excitation energy relative to other intermedi- 
ate transition states formed during the collision. Our 
choice of the reactions mentioned above was motivated 
by the fact that the former system has a positive Qgg 
value (16.6 MeV), whereas the latter system has a nega- 
tive one (-14.2 MeV). 

In Fig. [T] we show the ion-ion interaction potential 
V(R) for the head-on (zero impact parameter) collision of 
the ^^0-1-^^0 system at various center-of-mass energies. 
These results are essentially the same as those published 
in Ref. except for the energy dependence of V{R). 
This dependence arises from the time available for the 
system to undergo rearrangements and partial equilibra- 
tion, which approaches the frozen-density limit at high 
energies [13 • On the same figure we have also shown the 
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FIG. 2: Excitation energy for the head-on collision of the 
^^0-1-^^0 system for various -Ec.m. values (solid curves) as 
a function of the ion-ion distance R. Also, shown are the 
corresponding collective kinetic energy Ekin values (dashed 
curves) . 

relative location of the '^^S ground-state binding energy 
as represented by the —Qgg value. In Fig. [2] the cor- 
responding excitation energies are shown as a function 
of R calculated from Eq. ([T]). The fact that the excita- 
tion energy should be asymptotically zero is a good test 
of numerical accuracy, which is reproduced quite well as 
can been seen from the figure. The final value of E* and 
smallest R value are chosen to be the ones correspond- 
ing to the smallest relative velocity or smallest collective 
kinetic energy in the vicinity of the potential minimum, 
which alternately can be referred to as the capture point. 
Naturally, some of these quantities can only be pinned 
down within the limits of density-constraint frequency. 
For sake of completeness in Fig. [3] we also show the long- 
time behavior of the potential and the excitation energy 
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FIG. 3: Long-time evolution of the excitation energy, E* , and 
the ion- ion potential, V{R), for the head-on collision of the 
16q_,_16q gygtem at Ecm. = 11 MeV as a function of the 
ion-ion distance R. 



for the ^^O+^^O system at E^.m. = 11 MeV. The fig- 
ure demonstrates very nicely that the majority of the 
entrance energy is absorbed into intrinsic motion such 
that the compound stage is stuck in the vicinity of a 
certain R value with rather small oscillations. In the en- 
trance phase of the collision, the excitation energy E* in- 
creases monotonically with decreasing R. After the point 
of closest approach has been reached, the system is stuck 
close to that distance and the dynamical evolution turns 
abruptly to irregular oscillations in E* as seen in Fig. [31 
which the typical energy fluctuations of an excited en- 
semble. These fluctuations are here rather large due to 
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FIG. 4: Excitation energy for the head-on collision of the 

16q_,_16q gyg^gjjj foj. 

VciriOUS -C/c.m. values at the point of cap- 
ture (solid blue line) and the excitation energy calculated from 
Eq. ((5]) (solid red curve). Other curves are described in the 
manuscript. 



the small particle number 
the energy width is AE « y^lGsFE* /{Ntt'^ 
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6 MeV, 



which fits nicely to the observed fluctuations. 

Next we wanted to demonstrate the conjecture that the 
value of the excitation energy measured at the capture 
point in TDHF is the excitation relative to the compos- 
ite or dinuclear system formed during the collision. In 
Fig. m we plot the excitation energy as a function of the 
center-of-mass energy for both the analytic expression 
of Eq. (O (solid red line) and the TDHF results at the 
capture point as discussed above (solid blue line). As 
expected the TDHF result is below the one generated 
from Eq. The difference between the two curves is 
shown by the dotted line. We have also calculated the 
energy difference of the composite system relative to the 
ground state, which is simply the Ejjc at the point of 
capture minus the ground state energy of the "^^S system 
obtained by an unconstrained Hartree-Fock calculation. 
This is shown by the black dashed curve. If we shift 
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FIG. 5: Internuclear potential for the head-on collision of 
the ""Ca-h'^'Ca system for various Ec.m. values. The relative 
ground state binding energy of the ^''Zr nucleus is represented 
by the —Qgg value. 

the TDHF result by these differences we obtain the blue 
dashed curve, which is in agreement with the result ob- 
tained from Eq. ([S]). We should emphasize at this point 
that for many reactions the excitation energy at the cap- 
ture point is of great interest as opposed to the ground 
state to ground state value, since after the capture point 
many different reaction possibilities exist. 

We have repeated the above study for the ^*'Ca-|-''°Ca 
system, for which the Qgg value is -14.2 MeV. In Fig. [5] 
we show the ion-ion potentials calculated using the DC- 
TDHF method at a set of center-of-mass energies. The 
energy dependence is analogous to the ^^0-1-^^0 case. 
Fig. [6] shows the corresponding excitation energies cal- 
culated via TDHF using Eq. Again, the excitation 
energies gradually rise as the nuclei enter the interac- 
tion regime while the collective kinetic energies show a 
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rise when the nuclei first experience the nuclear attrac- 
tion but eventually fall due to slowdown of the composite 
system. Finally, we again plot the center-of-mass energy 
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FIG. 6: Excitation energy for the head-on collision of the 
'"'Ca-|-^''Ca system for various -Bc.m. values (solid curves) as 
a function of the ion-ion distance R. Also, shown are the 
corresponding collective kinetic energy Eum values (dashed 
curves) . 

dependence of the excitation energy at the point of cap- 
ture in Fig. [7] (blue solid curve) together with the one 
obtained from Eq. ([5]), except this time using the nega- 
tive Qgg value (solid red line). Again the two curves run 
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FIG. 7: Excitation energy for the head-on collision of the 
^^Ca-f^^Ca system for various Ec.m. values at the point of 
capture (solid blue line) and the excitation energy calculated 
from Eq. (O (solid red curve). Other curves are described in 
the manuscript. 

parallel to each other and their difference is shown by the 
dotted line. The long dashed line shows the difference in 
energy between the density constrained energy Euc 



the point of capture minus the ground state energy of the 
^'^Zr nucleus. The two curves are almost exactly the same 
and shifting the TDHF excitation curve by this energy 
difference produces the blue dashed curve, which is in 
excellent agreement with the one obtained using Eq. ([5]). 
Thus, we can safely say that the microscopically calcu- 
lated excitation energy represents the excitation energy 
at the point of capture. 
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FIG. 8; Internuclear potential for the head-on collision of the 
i6Q_|_34j,jg system for various i?c.m. values and two alignments 
of the '^'*Ne nucleus. The solid lines denote the potential for 
vertical alignment whereas the dashed curves are for the hor- 
izontal alignment of the "^^Ne nucleus with respect to the col- 
lision axis. 

It is also possible to obtain an approximate temper- 
ature by relating the excitation energy to the tempera- 
ture using the Fermi gas model, E* = aT^, where a « 
A/8 MeV~^ is the level density parameter. For the above 
reactions and the Ec.m. values used, this translates into 
a temperature range of 2.0 - 3.5 MeV for the ^^O-h^^O 
system and a temperature range of 1.4 — 2.5 MeV for the 
^"Ca-f^'^Ca system. The reliability of the above approx- 
imation should be higher for heavier systems. 

In order to examine how deformation influences the ex- 
citation energy during a heavy-ion collision we have cho- 
sen to investigate the ^^0-|-'^''Ne system. In Hartree-Fock 
calculations the neutron-rich '^'*Ne nucleus has a large ax- 
ially symmetric deformation. In the past we have exam- 
ined the effects of deformation on the ion-ion potentials 
due to the different initial alignments of the deformed 
nucleus [13, [2l|. We have performed TDHF collisions 
of ^^0-1- ^"^Ne at various center-of-mass energies and for 
the two extreme alignments of the "^^Ne nucleus, one in 
which the elongation axis is along the collision axis and 
the other for which it is perpendicular. In Fig. [5] we show 
the ion-ion potentials obtained for this system using the 
DC-TDHF method. We observe that at Ec,„,, = 11 MeV 
the vertical alignment case does not fuse, whereas the 
horizontal alignment does. Furthermore, the centroids of 
the potentials are different for the two alignments since 
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the nuclei start to come into contact earlier/later for hor- 
izontal/vertical alignments. Otherwise the energy depen- 
dence of the potentials is commensurate with the spher- 
ical systems. The fact that different orientations of the 
deformed nucleus lead to a difference in the time of con- 
tact between the target and the projectile is expected to 
manifest itself in the evolution of the excitation energy 
as well. Indeed, in Fig. [5] we see that the excitation for 
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FIG. 9: Excitation energy for the head-on collision of the 
i6Q_|_34-|yjg system for various Ec.m. values as a function of 
the ion-ion distance R. The solid lines denote the excitation 
energy for vertical alignment whereas the dashed curves are 
for the horizontal alignment of the ^*Ne nucleus with respect 
to the collision axis. 



the horizontal alignment of the Ne nucleus starts ear- 
lier but somehow does not reach as large a value as the 
vertical alignment case. This is a very interesting obser- 
vation since it would indicate that due to the differences 
in the excitation energy at the point of capture differ- 
ent alignments will likely have different probabilities for 
various exit channels. For example, this could be a very 
important consideration for superheavy formations. 

In summary, we have outlined a microscopic approach 
for calculating excitation energies of composite or din- 
uclear systems formed during heavy-ion collisions. The 
goal of the approach is to provide estimates for excita- 
tion energies at the initial point of capture, after which 
a multitude of exit-channel possibilities may exist and 
will be strongly influenced by the amount of excitation. 
The premise of our approach depends on the reliability 
of TDHF theory for describing the early stages of heavy- 
ion reactions. As discussed in the manuscript, there is 
mounting evidence that TDHF does provide a reliable 
description of heavy-ion dynamics and dissipation in this 
stage whereas the long-time evolution may be question- 
able. In order to elucidate the above arguments we have 
performed a number of calculations involving both spher- 
ical and deformed nuclei. It is our long-term goal to 
extend these calculations to heavier systems and in par- 
ticular to superheavy nuclei. 
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